November 30, 1994 HU Berlin-IEP-94/19 

Bielefeld BI-TP-94/47 



Polyakov loop and spin correlators on finite lattices q 

A study beyond the mass gap 

J. Engels 1 , V. K. Mitrjushkin 2 f\ and T. Neuhaus 1 

1 Fakultdt fur Physik, Universitdt Bielefeld, 33615 Bielefeld, Germany 

2 Fachbereich Physik, Humboldt- Universitdt, 10099 Berlin, Germany 

Abstract 

We derive an analytic expression for point-to-point correlation functions of 
the Polyakov loop based on the transfer matrix formalism. For the 2d Ising model 
we show that the results deduced from point-point spin correlators are coinciding 
with those from zero momentum correlators. We investigate the contributions 
from eigenvalues of the transfer matrix beyond the mass gap and discuss the 
limitations and possibilities of such an analysis. The finite size behaviour of 
the obtained 2d Ising model matrix elements is examined. The point-to-point 
correlator formula is then applied to Polyakov loop data in finite temperature 
SU(2) gauge theory. The leading matrix element shows all expected scaling 
properties. Just above the critical point we find a Debye screening mass Hd/T m 
4 , independent of the volume. 
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1 Introduction 



The determination of the correlation length £ and the screening mass /i£> from point- 
to-point correlation functions of the Polyakov loop is a non-trivial task, especially close 
to the critical point of lattice gauge theories. The difficulties are resulting on one hand 
from finite volume effects due to the nearby transition and on the other hand from the 
unknown parametrization of the heavy quark potential in the non-perturbative regime. 

In the transfer matrix (TM) formalism the levels of the transfer matrix provide an 
access to both £ and without the introduction of an ansatz for the quark potential. 
The formalism was first applied with success to the exact solution of the 2d Ising 
model [II], . It proved to be very efficient as well in the analysis of the 4d Ising model 
II, |J and also in the investigation of zero momentum correlators in Yang Mills theories 
P, Q . The zero momentum ( or plane-plane ) correlation functions are very convenient 
quantities to evaluate, since their TM form is simply exponential and the levels and 
matrix elements of the transfer matrix may be obtained easily from fits. 

To make the advantages of the TM formalism also available for the analysis of 
point-to-point correlation functions on d— dimensional spatial lattices we derive in this 
paper the corresponding TM expressions. In the 2d Ising model we test and confirm 
then the validity of our TM formula by comparison to the results obtained from plane- 
plane correlators. Simultaneously we are able to determine, where - as a function of 
the coupling constant - levels beyond the mass gap are of importance and what can be 
expected from such an analysis. 

We apply then the TM technique also to the case of the (3+1) dimensional SU{2) 
gauge theory. Our special point of interest is here, in contrast to the intention of 
an earlier paper |7j on the subject, the study of the higher levels, in particular their 
connection to the screening mass //d, and their influence on the determination of the 
mass gap. 

Let us consider d— dimensional spatial lattices with periodic boundary conditions 
of size iV d_1 L, where N denotes the number of points in each transverse direction and 
L that in one selected direction ( the z— direction ). The lattice spacing a is set to 
unity in the following. The partition funtion is then 



and V is the transfer matrix in z— direction. Its eigenstates | n) are chosen to be 
orthonormal. We order them such that we have for the eigenvalues A n (n = 0, 1, 2, ...) 



In addition we normalize our partition function so that we have for the vacuum state 




(1) 



V \n) = X n - \n) ; X n = e 



(2) 



Ho < Hi < fi 2 < ■■■ ■ 



(3) 



A = 1, /io = 0. 



(4) 
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This implies 

z = EH vL N 

n 

= 1 + e^ lL + e~^ L + e~^ L + . . . , 
and the partition function varies then between 1 and 2. 

2 Plane— plane and point— point correlators 

2.1 Plane— plane correlators 

We define zero momentum operators by 



(5) 



0(z) = N-^-Y,0(x ± ,z) , 



(6) 



XJ_ 



where 0(x±,z) is the Polyakov loop V(x±,z) for the 3 + 1 dimensional SU(2) 
gauge theory (d = 3) and the spin a XyZ for the two dimensional Ising model (d — 2). 
The corresponding correlation functions are 



F(z) = (O(z)-6(0)) = Z" 1 • Tr [(5(0) • V* • (5(0) • V 



z- 1 ■ E c 



.2 

mn 



e -VmZ . e -jx n (L-z) _|_ g-/J„z . e -jx m (L-z) 



n<m 



Z ~ K E 



„2 _ p -ii n L 



n<m 



where 



(n| 0(0) |m) , 



(7) 



are the level differences and the transition matrix elements. Because the eigenstates 
of the transfer matrix are either symmetric or antisymmetric under tranformations, 
which change the sign of O, we have c nn = . 

Let us compare the contributions of the different states to the correlator T(z) 
above and below the critical point f3 c , where /3 is the coupling (to be specified 
below). To do this we show explicitly the contributions to T(z) associated with the 
three lowest nonzero states 
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T(z) 



'10 



-fJiXZ 



+ e 



-Hi(L—z) 



+ C 



30 



g-A«3^ _|_ e -/i3{L-z) 



+ c 21 • e 



e -M21^ _|_ e -M2l(i-«) 



+ ... } 



(9) 



2.2 Point— point correlators 



In our subsequent analysis we establish a connection between the plane-plane and 
point-point correlators 

r(f) = (0(f) • o(o)) , 



(10) 



in the context of the transfer matrix formalism. 

The Fourier transform of the plane-plane correlator 



is, using eq.(|7D 



f(p z )=£e^.I>), 

z 

f(p z ) = Z- 1 .J2c 2 mn -e~^ L -G( Pz ^ mn ) 



where 



G(p z ; fj) = 2 (l - e~^ L ) sinh/i 



i -i 



4 sinh 2 - + 4 sin' 

2 2 

is just the Fourier transform of the sum of exponentials [e~^ z + e~^ L ~ z ^] 
connection between T(z) and the point-point correlator T(xj_;z) = T(x) is 

f(z)=J2T(x ± ;z) , 

x± 

and their Fourier transforms are related by 

f(p*) = T (f± = 0>P«) • 
To obtain the full correlator r(pj_,p z ) = T(p) we use the substitution 



A • 2 P z 

4 sin — 



4 sin - 



(11) 
(12) 

(13) 
The 

(14) 

(15) 
(16) 



where -D(p) is the lattice laplacian in <i-dimensional momentum space. Therefore we 
arrive at 



2 _ p ~^nL 



(17) 
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where 



G(p;n) =2(l-e~» L ) sinh/x 



4 sinh 2 - + D(p) 
2 



-i 



and 



G( Pz ;fi) = G(p ± = 0, Pz ;fi) . (19) 

Finally, to get the correlator T(x) we perform the inverse Fourier transformation, 
resulting in 



m = z ~ 1 -{ 4 • g{ 



+ C3 -G(f;^ 3 ) 

+ c 2 21 -e-^ L -G(x;fi 21 ) + ... } 
and G(x; /i) is just the Fourier transform of G(p; fi) 



G(x; n) 



1 



N*- X L 



<£ e -*-P. G{f]fl) . 



(20) 



(21) 



The ansatz we used in eq.(|16|) to obtain the point-point correlator is so far without 
proof. Yet the results of our numerical analysis of correlators strongly support it (see 
below) . 

An expression for the expectation value of the square of the lattice average of 
the operator 0(x) in terms of the matrix elements c mn may be easily derived in the 
following way. Since 



we obtain 



with 



N d -'L{o 2 ) =^r(f) = r(p = o), 

X 

N^LiO 2 ) =Z- 1 -J2 cln ■ e^ L ■ G{p= 0; ^ 



G(p = 0;n) = (l-e^ L ) coth 



:22) 



(23) 



(24) 



The correlator T(x) can be represented in the form of a superposition of two 
Yukawa- type potentials only in the case of the large- volume limit. If j\i-N ~ 1 then 
the finite-volume corrections are too strong and such a representation is not possible. 
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2.3 Correlation length and scaling behaviour 

Below the critical point (3 < (3 C the lowest nonzero energy level Hi - the mass gap 
- determines the large distance behaviour of the correlation functions. We therefore 
define the correlation length at (3 ~ (3 C by (see also ref.0) 

UP) = ^ ~\P-Pc\- V ■ (25) 

The contribution of the next level with energy /i 2 is assumed to be suppressed because 
of the additional factor e _MlL (or e~^ 2L ) in the third term of the right hand sides of 
eqs.(^) and Q2D|). Therefore the third level with energy /i 3 gives the next to leading 
corrections at large distances. 

The situation is different at (3 > (3 C . There the mass gap fix ~ , if iV is large 
enough and L ~ iV .In this case the first term on the right hand sides of eqs. (|]) and 
(P0|) becomes z — independent and the large distance behaviour is given by the next 
level difference A/i = /i 2 i ~ /! 2 or /i 30 = /i3 , so that the Debye mass is 

tiD = Afi. (26) 

Here //d is the nonperturbative equivalent to 2mo, where tub is the perturbative 
screening mass. Near the phase transition point /3 ~ (3 C all three levels are expected 
to give an essential contribution. 

In the thermodynamic limit below the phase transition point, j3 < /3 C , the correlator 
T(x) decays exponentially at large distances |x|3> 1 

r(aO~exp(- | f| /£-(/?)), (27) 

which entails 

Y,T(x ± ;z)<°o. (28) 

From eq. ( p8|) we conclude that for a finite lattice size all matrix elements are indepen- 
dent on iV in the large volume limit so that 

4 ~ N ° ! N - oo. (29) 

Well above the transition point, (3 > (3 C , the behaviour of the correlator T(x) for 
large separations \x\^ 1 is 

T(x) - a ~ exp(- \ x\ /Z+(p)) ; -> oo, (30) 

where a is a positive constant, independent of |x| and of for large but finite A^ 
so that 

J2r(x ± ;z) = f(z)~N d -\ (31) 

XJ_ 
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The constant is connected to the existence of a nonvanishing spontaneous magnetiza- 
tion and is in fact equal to (O 2 ) on finite lattices [[J . As the major contribution 
to T(z) is coming from the first term in eq.(^), which is proportional to cf and 
essentially independent of z , since \L\L <C 1 . We expect therefore that 



cf ~ N d - L ; N -f oo. (32) 

The same conclusion may be drawn as well directly by assuming A^— independence of 
{O 2 ) and eqs. (||) and (H). 



In the very neighbourhood (3 ~ f3 c of the phase transition we may apply finite size 
scaling techniques || [Uj to derive the A^— dependence of the matrix elements. Let 
us assume for simplicity, that N = L . According to finite size scaling theory any 
observable O with critical behaviour is supposed to have the following form 

O = N p ' u ■ foixN 1 ^) ■ N -f oo, (33) 

for fixed small x = ((3 — /3 C ) / f3 c . Here p is the critical exponent of the observable O 
and v the one of the correlation length. Due to eq. (p5|) , we expect then that the mass 
gap behaves for (3 < (3 C (x < 0) as 

fii = N" 1 ■ UixN 1 /") ; x ~ . (34) 
In the same (3— region the susceptibility x is defined as follows 

X = N d (0 2 ) ; (3 < p c . (35) 
Its critical behaviour is governd by the exponent 7 so that 

N d (0 2 ) = W tv ■ f x {xN x / v ) ; x ~ . (36) 

Since again the leading contribution to (O 2 ) is proportional to the matrix element 
cf Q we find after combining the last equation with eq.'s (|23|), (|24|) and ( p4[ ) that 

cj = N- 1+ ^ ■ UxN 1 '") ■ x ~ . (37) 

Although the latter equation was formally derived only for (3 < (3 C , we expect, because 
of the analytic (3— dependence on finite lattices, that it will be valid also above the 
critical point. 



3 Numerical results 

3.1 The twodimensional Ising model 

The 2d Ising model provides an ideal test case for a comparison of point-point and 
plane-plane correlators. This is so, because the levels \x n are all explicitly known [0. |2|]. 
Therefore a fit of the correlators in terms of the TM formulae, eqs. (||) and (|20D , requires 



7 



only the determination of the matrix elements. Moreover, Monte Carlo simulations of 
the model are relatively simple. 

To be more specific, consider a d — 2 Ising system on a lattice of size L x -L z = N-L 
with periodic boundary conditions. At every site i = (x, z) there is a spin a, = ±1 . 
The partition function Z is of the form 

Z= £ expjV £ OioA , (38) 

{<7fe=±l} V <V> I 

where f3 is the inverse temperature and < ij > means that only nearest neighbours 
interact . 

The eigenvectors | n) of the transfer matrix as well as the eigenvalues A n 
correspond to different numbers of collective excitations (quasiparticles) in transverse 
direction with momenta q%, q^-, ■ ■ ■ ■ These momenta take values (for even N) 

7T 3tt 7r(JV-l) r 

q = ±— , ± — , ± — for an even number of quasiparticles 

N N N 

2tt vr(iV-2) 
q = 0, ■■■5 i > n f° r an °dd number of quasiparticles. 

We may characterize the eigenvectors | n) via the momenta of the quasiparticles 
qi < q-i < . . . < q m as \ki,k 2 , . . . ,k m ) where q = j?k . Of course, the eigenvectors 
and eigenvalues depend on iV but do not depend on L . 

The four lowest states are then the vacuum | 0) , 1 1) = | k\ — 0) , | 2) = 
— |; |) and | 3) = | —1; 0; 1) . Above the vacuum level fi = , the next smallest 
jj, - the mass gap \L\ - is nonzero in the limit iV — > oo below the critical point 
(3 < fi c . At the critical point /xi(/? c ) ~ A^" 1 , as expected from eq . (|34]) , and above 
the critical point Hi tends to zero with increasing iV as ji\ ~ 

N -i/2 e -*-N _ In the 

thermodynamic limit the vacuum becomes therefore degenerate. In Fig. 1 we show for 
N = 30 the dependence on (3/ '/3 C of the levels /ii to /i 5 and the smallest and therefore 
most relevant level differences, which appear in the correlator formulae. 

To test the substitution, eq. (|l6|) , and the resulting eq.( |^) for the correlator T(x) , 
we have measured plane-plane and point-point correlators of the spin operator Oi on 
N = L = 30, 40, 50, 60 lattices. At each point 500000 cluster updates were performed 
and measurements taken every 10th update. Subsequently we have carried out fits to 
both correlators with varying numbers of levels to obtain the matrix elements. The 
results for the matrix elements on an iV = L = 30 lattice are compared in Fig. 2. In 
each case we show fits including either all distances r =| x |= 1, 2, 15 or only those 
for r > 2 . If more levels than assumed in the fit are contributing, we see a dependence 
on the lowest distance r taken into account in the fit. The effect is more pronounced 
for the point-point correlator than for the plane-plane correlator. Summarizing the 
experiences we have made with the different fits, we observe, that both formulae lead to 
exactly the same results, whenever the maximal number of levels is taken into account, 
which lead to non-negative (? mn , i.e. our ansatz is definitely confirmed. 
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The final result for the 2c? Ising model and N = L = 30 is shown in Fig. 3. We 
find that for (3 < 0.92/3 c only one term with /xio = /xi, the mass gap, contributes; near 
(3 C up to three terms are essential and well above the critical point, for (3 > 1.1/3 C , 
where the mass gap /xi ~ only one more term is present. The matrix element is 
increasing below (3 C with decreasing (3. However, the relevant factor in the correlator 
formulae, exp(— /XiL) , is negligible below (3 = 0.92/3 c , so that can no longer 
be determined from the fits. 

We have also studied the scaling properties of the major matrix element c{ in 
the three different (3— regions. As can be seen from Fig. 4, the predictions, eqs. fl29|) , 
(|3^) and (|37|) are all nicely confirmed by our results. Here we have used the known 
2c? Ising model values for 7 = 7/4 and v = 1. Obviously, we could not check the 
iV— independence according to eq.(p9[) of the higher matrix elements such as j 
since - as mentioned - their contributions are negligible for (3 well below (3 C . 



3.2 SU(2) lattice gauge theory 

We now want to apply our TM formula for the point-point correlator to SU (2) gauge 
theory. We consider a finite lattice of size L t ■ L d ^ x ■ L z = L t ■ N^ 1 ■ L {d = 3) with 
periodic boundary conditions. The standard Wilson action for SU{2) gauge theory is 

S w = P ■ E - \ TrU *) > ( 39 ) 

where (3 = 4/g 2 and □ = (x; up) refers to location and orientation of the plaquette. 
The field variables U K (x;t) £ SU(2) are defined on the links, and the U a are 
plaquette variables U □ = U X . K „ = U x . K U x+K . p Ul +p . K Ul . For the Wilson action the 



transfer matrix V is proven to be positive definite |TT|, [12[. Also, due to the Perron- 
Frobenius theorem [^] the vacuum | 0) is unique for a finite system, and we can 
choose again the normalization such that /xq = . 

Here, the Polyakov loop V(x) takes the role of the operator 0{x) . It is defined 
as usual by 



V(x) = ^Tr 



Lt 

]ju A (x } t) 

U=l 



(40) 

The Monte Carlo data for the point-point correlator, which we want to analyze 



in the following, were computed |14| on L t = 4, iV = L = 12, 18, 26 lattices with 
10 5 — 4 • 10 5 updates and measurements every 10th sweep. In contrast to the case of 
the 2d Ising model, in SU (2) gauge theory the level differences are unknown and have, 
like the matrix elements, to be determined through the fit. 

In general we find a behaviour resembling very much the one of the 2d Ising model. 
In particular, the number of levels, which may be extracted from the fits is comparable. 
Fits with more than two levels are only possible on the largest lattice very close to the 
transition. Otherwise one either obtains negative squares of matrix elements or there 
is no minimum of x 2 ■ Taking into account more than one term in eq. (|20|) tends to 
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decrease the result for the mass gap level This can be seen in Fig. 5, where we 
show the results for N \x\ from one and two level fits. The inclusion of a third level in 
the fit, however, does not change \x\ anymore. 

The fit result cf Q /Z for the major matrix element was subsequently checked for its 
scaling properties. Note, that here the partition function Z is not explicitly calculable 
from eq. ([5]), because we know only the lowest level(s). Like in the case of the 2d Ising 
model we find all predictions from eqs. (p9|), (|3~2"D and (|3~7D very well confirmed. This 
is shown in Fig. 6. In the scaling test, Fig. 6b, we have used 7 = 1.24 and v = 0.63, 



the values of the 3c? Ising model, in accord with the universality hypothesis |15| , which 
predicts equal critical exponents for SU(2) and the 3d Ising model. 

It is interesting to look at the behaviour of the next to leading level ( or level 
difference ) A^u. As can be seen from Fig. 7, Afi drops from a higher value below f3 c 
at the transition to a value near to one ( in lattice units ) and stays then relatively 
constant and moreover independent of the lattice sizes used here. This second level 
fixes the large distance behaviour of the correlation functions above /3 C , since fix, as 
is evident from Fig. 5, is essentially zero there and a third level does not contribute 
outside the transition region. Therefore we identify it with /i£>. Because we have L t = 4 
and T = 1/L t we are led to a ratio fio/T ~ 4, slightly higher than the ratio found with 
conventional methods[16|]. On the other hand a higher value seems to be preferred by 



next-to-leading order perturbation theory calculations [17 . In the close vicinity of the 



transition we have found at four j3— values solutions to three level fits on the N = 26 
lattice. They are also shown in Fig. 7. We see that the lower of the two levels beyond 
the mass gap is approaching zero at the transition. Indeed, if interpreted as fin, the 
expected iV— behaviour at the transition is propotional to N^ 1 . 

Finally we present in Fig. 8 the second matrix element c\/Z resulting from two 
level fits. The fluctuations of c^/Z in the neighbourhood of the transition are proba- 
bly due to both the statistical uncertainties in the data and the possible influence of 
higher levels. It is remarkable, that outside the transition region the matrix element is 
iV— independent, constant and moreover about equal well above and below the critical 
point. Compared to the first matrix element, however, the second one is rather small. 



4 Summary 

We have derived a transfer matrix formula for the point-to-point correlation function 
on d— dimensional spatial lattices. The advantage of such an approach lies in its direct 
access to the correlation length and/or the screening mass. Moreover, the disconnected 
point-point correlation functions may be analysed without the need for any subtraction. 

The formula was tested in the 2d Ising model by comparison with plane-plane 
correlation functions. We find that from both observables the same information may 
be extracted, whenever the maximal number of levels are taken into account, which lead 
to physically meaningful results in the fitting procedure. Quite naturally a difference 
is observed, if existing higher level contributions are neglected. 
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In both the 2d Ising model and the (3 + 1) dimensional SU(2) gauge theory we 
found the same general behaviour of the levels and matrix elements : 

1. well below the critical point only the mass gap fii is contributing to the correla- 
tors, the matrix element cf is independent of N ; 

2. close to the transition up to three levels are contributing and cf ~ N^ 1 +f/ u • 

3. well above the critical point only one higher level beyond the essentially zero 
mass gap is contributing and c^g ~ 

N d-i _ 

The detection of still higher levels seems to require an extremely large statistics of 
the data and very large lattice volumes. Most probably there is only a chance for such 
a program close to the transition point. 

In SU (2) gauge theory we have calculated the change in the mass gap due to the 
presence of the higher levels. We find that this effect decreases the mass gap result. 
Finally we have determined the Debye screening mass to yU D /T 4, independent of 
the lattice size used. 
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Figure captions. 



Fig.l The lowest energy levels in the 2d Ising model as a function of (3//3 c f° r N = 30. The 
mass gap fxi and (12 up to ^5 are shown as solid lines, the differences ^21 = fJ-2 — Mi and 
M32 = M3 — M2 as dashed lines. 

Fig. 2 Comparison of matrix element fits with one level (a), two levels (b) and three levels 
(c) in the N = L = 30 2d Ising model vs. P/f3 c . The results shown are from fits to zero 
momentum correlators including all distances for r > (dashed lines) and r > 2 (solid 
lines). The corresponding results for point-to-point correlators are shown by dotted and 
dashed-dotted lines. 

Fig.3 The lowest levels fii, ^3 and M21 (dashed lines) and the corresponding best fit matrix 
elements Ci ,c^ and c^e - ^ as a function of /3/(3 c in the N = L = 30 2d Ising model. 
The dotted line is c^. 

Fig.4 Scaling properties of the mass gap matrix element cf below (3 C ( ~ iV° ), 
close to p c ( ~ AH^-i 

as a function of xN , x = ((3 — (3 C )/(3 C ) and above j3 c 
( ~ AT ) for N = L = 30,40,50,60 (solid,dashed,dashed-dotted,dotted lines) in the 2d 
Ising model. 

Fig.5 The dependence of iV/xi on the number of levels used in the fit for N = 18 (squares) 
and iV = 26 (diamonds) as a function of /? in SU{2) gauge theory. Two level fits are 
shown by filled symbols, one level fits by empty ones. The inset shows for N = 26 also three 
level fits (circles). 

Fig. 6 Scaling properties of the mass gap matrix element c±q/Z in SU{2) gauge theory. 
The figure corresponds to Fig. 4. Here N = 12, 18, 26 (crosses, squares,diamonds). 

Fig.7 The level difference A/j, from two level fits in SU(2) gauge theory as a function 
of P . The notation is the same as in Fig. 6. The results of three level fits on the N = 26 
lattice for the two levels beyond \i\ are shown as filled diamonds. 

Fig.8 The second matrix element c^jZ from two level fits in SU{2) gauge theory as a 
function of (3 . The notation is the same as in Fig. 6. 



13 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94 1 2003 v 1 



2 . 



1 . 5 



. 5 



. 





I I I 

N = 30 


A \ 


- 


\ u \ \ 
\ 4 \ \ 






^ / /A/ 


\ \ 

N \ V 

S N \ \ 




V x \ \ 

NS ^ \ \ 
X- x \ \ 




- >%- . . \ v 

\> \ \ 
\\ ^ \ 




\. v. -- _\ 


^32 


I I I 


i i i 



0.7 0.8 0.9 1.0 1.1 1.2 1.3 



Figure 1 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94 1 2003 v 1 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94 1 2003 v 1 




Figure 3 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94 1 2003 v 1 



This figure "figl-5.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/94 1 2003 v 1 




2 . 25 



2 . 30 



2 . 35 



Figure 5 
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Figure 7 
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Figure 8 



